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ABSTRACT 

White dwarfs experience a thermal renaissance when they receive mass from a stellar companion in a binary. 
For accretion rates < 1O~ 8 M yr" 1 , the freshly accumulated hydrogen/helium envelope ignites in a thermally 
unstable manner that results in a classical novae (CN) outburst and ejection of material. We have undertaken a 
theoretical study of the impact of the accumulating envelope on the thermal state of the underlying white dwarf 
(WD). This has allowed us to find the equilibrium WD core temperatures (r e ), the classical nova ignition masses 
(Mgn) and the thermal luminosities for WDs accreting at rates of 10~ u - 1O~ 8 M yr" 1 . These accretion rates 
are most appropriate to WDs in cataclysmic variables (CVs) of P 01 -b < 7 hr, many of which accrete sporadically 
as dwarf novae. We have included 3 He in the accreted material at levels appropriate for CVs and find that it 
significantly modifies the CN ignition mass. We compare our results with several others from the CN literature 
and find that the inclusion of 3 He leads to lower M ign for (M) > 1O" 1O M yr" 1 , and that for (M) below this the 
particular author's assumption concerning T c , which we calculate consistently, is a determining factor. Initial 
comparisons of our CN ignition masses with measured ejected masses find reasonable agreement and point to 
ejection of material comparable to that accreted. 

Subject headings: binaries: close — novae, cataclysmic variables- nuclear reactions, nucleosynthesis, abun- 
dances — stars: dwarf novae — white dwarfs 



O 1. INTRODUCTION 

Classical novae (CN) are the result of unstable thermonu- 
5_i ■ clear ignition of freshly accreted hydrogen and helium on a 
\ white dwarf (WD) in a mass transferring binary. The orbital 
Ci ■ periods of these binaries presently range from 1.4 to more 
ILJ ' than 16 hours, with a large number in the 2-4 hour range 
. ' JDiaz & Bruchill997tlwirneill2002ft . These cataclysmic vari- 
ables fCVs: lWarnerJI 1995ft have low-mass stellar companions 
and spend most of their lives accreting at time averaged rates, 
(M) w 10" 9 - lO" n M yr" 1 (e.g. Howell, Nelson, & Rappa- 
port 2001). The frequency of CN as a function of orbital pe- 
riod depends on the accumulated mass immediately prior to 
the explosion, M lgn . However, the value of M; gn strongly de- 
pends on the WD's core temperature, T c , which is a priori 
unknown and must either be calculated or measured. Knowl- 
edge of the ignition mass is also needed for comparisons to 
measured ejected masses in CN that would decide whether 
or not the underlying carbon/oxygen WD is being excavated 
during the CN. 

CVs are formed when the WD made during a common en- 
velope event finally comes into contact with its companion 
as a result of gravitational radiation losses over a few Gyr 
(e.g. Howell et al. 2001). The WD will have cooled during 
this ti me; a Q.2QM G) He WD would have T c = 3.3 x 10 6 K at 
4 Gyr ilAlthaus & Benvenutoll 19971) . whereas a 0.6M^ C/O 
WD would have T c = 2.5 x 10 6 K in 4 Gyr (Salaris et al. 
2000). These WDs would have effective temperatures T e ff ps 
4500-5000 K. A subset of the CVs, called dwarf novae 
(DN), contain a WD accreting at low time-averaged rates 
(M) < 1O~ 9 M yr" 1 , where the accretion disk is subject to 



a thermal instability which causes it to rapidly transfer mat- 
ter onto the WD (at M S> (M)) for roughly a week once ev- 
ery month to year. The M onto the WD is low enough be- 
tween outbursts that the UV emission is dominated by the in- 
ternal luminosit y of the WD , allowing for a measurement of 
r eff > 10,000 K lSionl ( 119991) . This is the best evidence that 
the WD is hotter than a coolin g WD of sim ilar age due to 
the thermal impact of accretion. ISionl (11995ft has argued that 
most of this heating is from gravitational energy released (i.e. 
converted to thermal energy) in the WD interior via material 
compression. This energy is then transported outward to the 
stellar surface, and we find that it can similarly heat the core. 
We will discuss this in great detail, as well as explain the im- 
portant role of nuclear "simmering". A thorough comparison 
of our theoretical work to T e g observations and a summary of 
the imp lications for CV evo lution are in a companion paper 
dTownslev & Bildstenl2003l) . 

The gravitational energy released (GM/R) when a particle 
falls from a large distance to the stellar surface, sometimes 
termed the accretion energy, is deposited at, or near, the pho- 
tosphere and is rapidly radiated away. This energy does not 
get taken into the star, as the time it takes the fluid to move 
inward is always much longer than the time it takes for heat 
to escape. This basically eliminates the outer boundary con- 
dition and instead points to the importance of energy released 
deep in the WD due to both gravitational energy released, via 
compression, within the star and hydrostatic nuclear burning. 
That energy takes a long time to exit, will still be visible when 
accretion has halted, sets T e ff, and reheats the WD. 

We begin in $2\by explaining the physics of heating within 
the accreted envelope and how that differs from heating in the 
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deep core (discussed in Appendix A). Since these WDs are ac- 
creting at rates low enough so that the accumulated hydrogen 
burns unstably, we also need to carefully track the accumu- 
lated mass of material and account for the additional energy 
source of slow nuclear "simmering" near the base of the ac- 
creted layer. This is significant as the layer accumulates and 
eventually triggers the CN. The calculation of M[ gn for a fixed 
T c is explained in ^3] (closely comparing to the original Fu- 
jimoto 1982 work), where we also exhibit the importance of 
the 3 He content in the accreting material. 

Our method of finding the equilibrium WD core tempera- 
ture (Townsley & Bildste n 2002) is explained in j4|and im- 
proves on that of Iben. Fuiimoto. & M acDonaldl ( 1992a) by 
allowing the accreted envelope mass to change through the 
CN cycle. We show that early in the cycle, when the mass 
of the newly accreted layer is small, compressional heating 
is small, and the WD cools. Later in the cycle, the accreted 
layer becomes thick enough that compressional heating along 
with nuclear simmering heats the core. So, the core cools at 
low accumulated masses and is heated prior to unstable ig- 
nition. Using CN ignition conditions to determine the max- 
imum mass of the overlying freshly accreted shell, we then 
find the steady-state (i.e. cooling equals heating throughout 
the CN cycle) core temperature, T cfiq , as a function of (M) 
and M. We close 50]by showing that for (M) < 10" 9 M Q yr" 1 
the WD reaches this equilibrium after accreting an amount of 
mass less than its own and that in CVs with P or b < 2 hours 
the accretion rate changes slowly enough that the WD should 
always be near the equilibrium given by (M) and M. 

Having found the core temperatures, we then predict the 
amount of mass that needs to accumulate on the WD in or- 
der to unstably ignite the accumulated hydrogen as a function 
of (M) and M. Our work is compared to previous calcula- 
tions in ^5] finding very large differences with others due to 
their lack of knowledge of T c . We also compare our predic- 
tions of Mj gn to those few measures of the ejected masses in 
CN of known orbital period. Within the confines of this lim- 
ited comparison, we find that the ejected masses are compa- 
rable to the amount accumulated. The equilibrium also yields 
powerful relations between the WD surface temperature and 
its accretion rate and mass. These allow for tests of binary 
evolution models as well as measures of the masses of these 
WDs; a comparison to the accumulated observations of effec- 
tive temperatures in the DN system s is being published sepa- 
rately (Townslev & Bildsten 2003). We close with a summary 
and discussion of future work, especially the seismology of 
these accreting WDs. 

2. ATMOSPHERIC MODEL INCLUDING ACCRETION 

Starting with Mestgjl (119521) single WD atmospheres are 
characterized by a relation between th e surface luminosity, 
L, and the core temperature, T r (e.g. IChabrier et a"fl 120001 
lHansenlU999t iSalaris et alJl20"Ool Moodily! " This is one 
of two components (the other being the thermal content of the 
core) needed for evolutionary calculations of a cooling WD. 
In this case the atmosphere extends to the location where elec- 
tron conduction forces the core to be isothermal. This occurs 
at M env = 5 x 1O _5 M for T c = 4 x 10 6 K and M = 0.6M Q when 
L = 3 x 10~ 4 Lq (see Figure 1 in Fontaine, Brassard, & Berg- 
eron 2001). 

In an accreting WD, energy is released in the atmosphere 
by the slow downward motion of matter in the gravitational 
field and by a slow rate of hydrogen fusion. This internal 



energy release modifies the L-T c relation, giving an accret- 
ing WD a higher L than a non-accretor of the same T c . The 
L-T c relation of non-accreting WDs depends on the surface 
gravity, g, which is set by the mass, M, and on the composi- 
tion of the envelope (which fixes the opacity), typically split 
into layers of pure H and He. The WD mass enters in much 
the same way in the accreting case, but the composition also 
determines the local energy release in the envelope. We pa- 
rameterize composition by the mass of the freshly accreted 
layer, M acc , which has the composition of the donor star. We 
use a solar composition with mass fractions: Xu = 0.709, 
X4 He = 0.271, X c = 0.009, and X = 0.011. Since 3 He will 
prove important as an ignition trigger (see discussion in ^3) at 
low (M), we will show models with Xi Ue = 0.001 and 0.005, 
the lower and u pper limits on expected values for the accreted 
mater ial (e.g. ID Antona & MazzitellHl 19821: llben & Tutukovl 
1984). Interior to this layer we assume a C/O core. As mate- 
rial is accreted at a time averaged rate (M), M acc varies from 
zero (immediately after a classical nova) to some maximum, 
Mi gn , just before ignition of the next classical nova. Thus three 
parameters, M, M acc and (M) lead to an L-T c relation. 

Due to the deposition of energy in the envelope, the sur- 
face luminosity is not equal to the thermal energy lost by the 
degenerate core. However, the atmospheric model provides 
not only L at the surface, but also the luminosity, L core , which 
enters the bottom of the atmosphere from the core. This lu- 
minosity couples the envelope model to the thermal content 
of the core, determining the temporal evolution of T c . As we 
will show, when T c is low, the envelope heats the degenerate 
core and L core < 0. 

2. 1 . Thin Shell Approximation 

As a spherical, gravitationally self-bound object accretes 
matter, it undergoes structural changes which we divide into 
two types: (1) global structural changes which affec t the core 
radius , and (2) internal movements of material. See lNomotol 
( 1982) and Appendix A for a discussion of the distinction be- 
tween these terms. In the outer parts of the WD, the two ef- 
fects are largely decoupled; material moves through the enve- 
lope and down into a deeper reservoir and the WD radius and 
the gravitional field in which the envelope exists are set by the 
total WD mass. This is particularly useful for the case of low 
(M), where the global changes happen extremely slowly, and 
can therefore be safely neglected. Our treatment of the core 
is presented in Appendix A, where we show that the energy 
release in the envelope dominates that in the core by a factor 
~5. 

Within a thin shell near the surface, accreted material enters 
at the top, and displaces pre-existing material which moves 
'out' of the shell, deeper into the star. With the assumption 
that the radius and total stellar mass, and thus g, are fixed 
on the timescale being considered, the local heat equation be- 
comes 

Ds ds ds dL 

T — =T— + Tv r —= + e N , (1) 

Dt dt dr dM r ' 

where is the nuclear burning rate, s is the entropy per gram, 
and v r = — (M) /4irr 2 p is the downward speed of material mo- 
tion in response to accretion. If a fluid element loses entropy 
as it is compressed, then that liberated energy is carried away 
by heat transport. The thermal structure of the envelope deter- 
mines the two spatial derivatives and the nuclear energy input 
in equation Q. Thus ds/dt, the time evolution of the en- 
tropy field, depends on the current thermal structure and (M) . 
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A thermal structure for which ds/dt = is a static state for 
the envelope. If a higher (M) were applied without chang- 
ing the thermal structure, a local buildup of heat would oc- 
cur, ds/dt > 0, due to the increase of \v r \ without a change in 
dL/dM r . A lower (M) leads to a corresponding cooling of the 
envelope. When the entropy advection term dominates over 
the nuclear burning term, this static state is stable because L 
has a stronger temperature dependence than s in both degener- 
ate and nondegenerate regimes. However, when e N dominates 
over the advection a thermal instability is likely, leading to a 
classical nova explosion. 

The radial derivative is evaluated from the envelope struc- 
ture. We work in pressure coordinates, and define m = 
(M) /A-Kr 1 , so 



ds . ds gmcpT 

Tv r — = SmT — = 

dr 6 dP P 



dlnT 



d\nP 



•V. 



ad 



(2) 



where cp = T{ds/dT)p is the specific heat at constant pressure 
and V a d = (dlnT /dlnP) s . The entropy advection is given by 
how much the temperature gradient differs from the adiabat. 
In the static state, the gradient of the luminosity itself is then 



dh 
dP 



= (M)c P 



dT 
dP 



dT 
~dP 



4irr 2 e N 



(3) 



which we use with the equations for heat transport and hydro- 
static equilibrium to solve for both the structure and luminos- 
ity in the envelope, resulting in an L and L C0le which depend 
on T c , M, (M), and M acc . The contribution to L from accretion 
is the heat not advected with the material, and is referred to in 
the literature as "compressional heating." Though this choice 
of words is unfortunate since the action occurring is more like 
cooling, we use this common term in our discussions. In all 
the calculations presented here we use a fixed R for the last 
term in equation 0, that of a zero temperature WD, 8.8, 5.7 
and 4.1 x 10 8 cm for M = 0.6, 1.0 and 1.2M Q respectively. 
As a single point check, matching the bottom of the envelope 
and the outer edge of a finite-temperature core model yeilds 
R = 8.9 x 10 8 cm for M = O.6M , (M) = 6 X 10- n M o vr" 1 , 
and M acc = 1.4 x 10" 4 M Q . By dropping the nuclear burning 
term and making a few other simplifying assumptions equa- 
tion can be integrated analytically to obtain an approxi- 
mate luminosity, L « 3(M)kT c / \im p (see AppendixlBl. where 
k is Boltzmann's constant, m p is the proton mass and p ~ 0.6 
is the mean molecular weight of the freshly accreted envelope. 

Since M acc changes with time as M acc (f) = (M) At, the static 
envelope structure is an approximation to the evolution. For 
this method to be accurate, we want the neglected term to be 
small, namely \ds/dt\ <C \v r ds/dr\. Evaluating at a fixed 
pressure and using equation (0 gives 



dT 
~dt 



dlnT 



dlnP 



(4) 



where r acc =P / grh is the accretion time. Dropping the temper- 
ature gradient (good in the conductive regions), finite differ- 
encing and setting At = r acc leads to the condition AT /T <C 
2/5, where AT is the temperature change at the envelope base 
during accumulation. As will be shown in Section|4] the evo- 
lution of M acc takes place at constant T c , so that this condition 
is always satisfied and our static structures are an excellent 
approximation for the accumulating envelope. 

Now we wish to consider the effect of the outer boundary 
condition on the structure of the envelope. To determine this 



we compare the thermal cooling time of an outer layer, t± = 
AMcpT/L to the time to accrete that layer r acc = AM/(M). 
Here AM is the mass of the (thin) outer layer being consid- 
ered, T is the temperature of the layer, and L is the overall 
luminosity which is driven mostly by energy released from 
deeper layers, L m c p (M)T c . Thus the ratio T t h/r acc mT/T c , 
so that where T <C T c in the outer layers, the material can cool 
much more quickly than it is being accreted and therefore has 
no influence over the thermal structure of deeper layers. In 
other words, the entropy of the outer envelope is determined 
by the luminosity flowing through it, rather than any exter- 
nal boundary condition. The condition T <C T c is sometimes 
violated during a dwarf nova e (DN) even t, where either the 
high external radiation field (Prinale 1988) or the more rapid 
accretion (Godon & Sion 2002i ISionl 1 9951) temporarily mod- 
ifies the outermost layers of the WD envelope. However, this 
heat does not penetrate to depths critical for the WD core, but 
only to a depth where the thermal time is of order the duration 
of the DN outburst, and is the likely cause of the fading T e s 
seen following a DN outburst (see Sion's 1999 review). 

In addition to the downward motion due to accretion, 
the heavier Helium ions will sink with respect to the over- 
all flow. While we have not included this effect, we esti- 
mate here when it is important. The drift velocity of the 
He ions is Vdrift — 2m n gD/kT where D is the th ermal dif- 
fusion coefficient (See Delove & Bildsten 2002). This is 
to be compared with the downward velocity of the fluid 
due to accretion, v acc = (M) /4TrR 2 p. In the simplest ap- 
proximation, Coulomb scattering off protons dominates and 



D = 0.024 A" 1 T, 



5/2 _! 

Pa 



where Tj = T/10 7 K and p4 
Coulomb logarithm. For Vdrift < v acc this yields the relation 



jAlcock&Tll a rionovl[T980h. 

p/10 4 g cm" 3 and A is the 



ship (M) > lO" n M yr _1 r 7 
-li 



3/2 



(M/0.6M Q )A" 
i 



so that for our 

lowest accretion rate, 1O~ U M yr" 1 , where Tj = 0.4, only the 
high mass models deserve further scrutiny. Below this ac- 
cretion rate, the energy release due to sedimentation of the 
Helium in the accreted material can be a moderate fraction 
of that released due to compression. For our M = 1 .0M Q 
model with (M) = 10 _11 Mq yr -1 , the conditions at the base 
of the H/He layer when M acc = 0.90M ign are p = 1.4 x 10 4 g 
cm" 3 , T = 4.1 x 10 6 K so that the Coulomb coupling param- 
eter r = 2. Diffusion coefficients in this gas-liquid crossover 
regime have yet to be authoritatively ca lculated. Using the 
fitting form of Iben & MacDonald ( 1985), which is very sim- 
ilar to the canonical results of iPaauet te et alj Q 9 86) in this 
regime (See Iben, Fuiimoto, & MacDonald 1992b), we obtain 
Vdrift/ Vacc = 1.3. However, using the liquid diffusion results of 
Wallenborn & Baus ( 1978), we obtain Vdrift/ V acc = 0.38. Re- 
gardless of this uncertainty, it is clear that in future work dif- 
fusion needs to be considered for M > 1 .0M Q and (M) < 
1 0" 1 ! M yr" 1 . As inferred from the amount of s edimentation 
found by Iben, Fuii moto. & MacDonald! ( fl99 2b). this limit is 
at higher (M) when a higher T c is assumed. 

2.2. Accreting Envelope 

We now present results from the envelope models, inves- 
tigating how (M) and M acc determine the L-T c relationship. 
The outer boundary is at P = 10 10 erg cm" 3 , in the upper edge 
of the radiative layers. We set the temperature there to that 
of a radiative zero solution for L = AirR 2 a^T^ f{ . Then equa- 
tion along with the other structure equations are integrated 
inward, keeping track of L(P) and T(P), to P = gM dcc /4irR 2 , 
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the base of the accreted layer. Interior to this, L = L core is 
taken to be constant and the integration is continued, using 
a pure carbon composition for simplicity, until P = 40Pi en , 
where Pi gn = gMi gn /AnR 2 and M[„ n is found consistently at 
each (M). A value of P = 2 x 10 20 erg cm" 3 is used for the 
M = O.6M demonstrations in this section, corresponding to 
the Mig n for an equilibrated model at (M) = 1O" 1O M yr" 1 . We 
label the temperature at the point where we stop the integra- 
tion T c , interior to which the core is isothermal. The equation 
of state for the ele ctrons and i ons w as taken from analytical 
approximations of Paczvhski ( 1983) with Coulomb correc- 
tions from Farouki & Hamaguchi ( 1993). The radiative opac- 
ities from OPAL ( Iglesias & Rogers 1996) and the conductivi- 
ties ofjltoh et al. ( 1983) were used. Nuclear re action rates and 
energies are from Caughlan & Fowler ( 1988) with screening 
enhancement factors for high density from lIchimaruN 19931) . 
which also appear in Ogata. Ivetomi. & Ichimaru ( 1991) 

The above treatment of the outer edge of the C/O core pro- 
vides a reasonable approximation to the evolution of this re- 
gion during a CN cycle without the complications of explicit 
simulation. At 40Pi gn the temperature change due to adiabatic 
compression is 1%. Under this approximation in the equi- 
librium state discussed in |4] fluid elements near the base of 
the accreted layer move on paths for which they must lose 
entropy during the CN cycle, i.e. ones more shallow than 
the adiabat in the T-P plane. Deeper fluid elements move on 
paths for which they must gain entropy, steeper than the adi- 
abat. This exchange is consistent with our expectation for the 
actual evolution, and since it only contributes peripherally by 
setting the gross temperature excursion from T c at the base of 
the accreted layer, the accuracy of a more sophisticated treat- 
ment is not considered worth its price in complexity at this 
point. The fraction of the core affected by the CN cycle in 
this treatment is also roughly consistent with that expected on 
energetic grounds. For example, a luminosity of 1O~ 3 L de- 
livers enough energy in the interval between classical novae, 
~ 10 6 years, to raise the temperature of the outer 1O~ 3 M of 
the core by 10 6 K (~ 10%), this mass coordinate corresponds 
to a pressure of just under 10 20 erg cm" 3 . 

The relationship between L and T c for several (M)'s are 
shown in Figure Q] for a O.6M WD (R = 8.76 x 10 8 cm) 
with M acc = 5 x 1O" 5 M . This M acc is less than the CN ig- 
nition mass for most of the range of (M) and is compara- 
ble to the hydrogen layers on cooling WDs. All the demon- 
stration curves in this section are for Xi He = 0.001. Shown 
for reference (dotted line) is the L-T c relation for an isolated 
O.6M WD with pu re layers of ma ss M H = 6 x 1O" 5 M and 
M He = 6 x 1O" 3 M (Hansen 1999). The corresponding re- 
lation for our envelope model, with (M) = 0, is shown as 
the dashed line and provides an initial comparison to previ- 
ous work. The offset from the cooling WD is due to the 
H/He envelope being more opaque than that of pure hydro- 
gen. The slightly larger divergence for 10 7 K < T c < 2 x 10 7 
K is likely due to the radiative-conductive transition point ly- 
ing below the outer H-rich layer, a location where the cool- 
ing model has a thick layer of helium whereas our enve- 
lope is pure carbon. The sharp change in the dotted line at 
L < 1O~ 4 L is due to surface convection extending down into 
the co nductive regions, commonly termed " convective cou- 
pling" (Fontaine, Brassard, & Bergeron 2001). To avoid this 
complication, we have restricted ourselves to L > 1O~ 4 L , 
where the outer convective layer need not be modelled di- 
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FIG. 1 . — Exiting surface luminosity, L, as a function of core temperature, 
T c , for a 0.6Mq WD. The accumulated H/He layer has a mass M acc = 5 x 
1O~ 5 M0. Various accretion rates are shown as labelled solid lines with (M) = 
10-", KT 10 - 4 , 1CT 10 , and 1O- 9 4 M yr 1 . The dashed line has (M) = 0, and 
the dotted line is the L-T c relation of an isolated cooling DA WD with pure 
layers of mass Mn = 6 X l(r 5 M Q and M Hc = 6 X 1(T 3 Mq |H ansen 1999). 



rectly to determine the L-T c relation. 

2.3. Parameter Dependences of White Dwarf Luminosities 

We now discuss how the L-T c relation changes under active 
accretion. It is only for low (M)'s and high T c 's that the lumi- 
nosity flowing through the envelope is dominated by the hot 
core. The T c below which the energy release due to accre- 
tion becomes important is determined by (M) ; higher values 
release energy faster via compressional heating and therefore 
their luminosity becomes dominated by accretion at a higher 
T c . Thus in Figure^ while there is a large area of the curve 
at (M) = 1O _11 M yr" 1 for which core cooling dominates at 
this M acc , there is little such area at (M) = 1O~ 94 M yr" 1 . The 
upturn in the luminosity at T c > 3 x 10 7 K is due to energy 
released from fusion of 3 He, a species not present in abun- 
dance in the cooling models. The onset of the nuclear energy 
source is quite sensitive to M acc since the conditions (T, p) at 
the base of the H/He layer dominate this effect (see Section|3] 
for details). 

The other major feature of the L-T c curve is the minimum 
L for a given (M) . In this situation the envelope determines 
its own temperature and the surface luminosity becomes in- 
dependent of T c . There is a temperature inversion in the en- 
velope: the temperature rises inward through the radiative 
layer, flattens out when conduction becomes dominant, passes 
through a T max and then decreases towards the core. The en- 
ergy liberated by compression goes into the core as well as 
coming out the surface. This lower limit, L(T C — * 0), is shown 
in Figure|2]as a function of (M). The upper solid curve is for 
M = O.6M and M acc = 5 x 1O" 5 M . Raising M acc does not 
change the relation - a curve for M acc = lO^M© is indistin- 
guishable - except to reduce the maximum L at which such a 
self-supported envelope state can exist. A higher M acc leads 
to a higher density at the base of the accreted layer and causes 
fusion to become important at lower envelope temperatures 
and luminosities. 
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FIG. 2. — Surface luminosity as a function of (M) for very cold WD 
cores. The upper solid line is for M = 0.6Mq and the lower solid line is for 
M = I.OMq both with M acc = 5 X 1CT 5 Mq. The dashed line is the values 
for the top solid line scaled by (1.0/0.6)~ 2 / 3 in order to demonstrate that 
our analytic WD mass dependence is close to correct. The dotted line is the 
analytical approximation discussed in the text. 
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FIG . 3 . — The WD surface luminosity L as a function of T c for various accu- 
mulated layer masses. The upper set of solid lines all have (M) = 10~'°' 4 Mq 
yr~' and have M acc = 0.5, 1, 3, and 5 X \(T^Mq as indicated. The lower set 
of solid lines has (M) = KT U M yr 1 and M acc = 3, 5, and 9 X 10" 4 */©. 
The WD mass is 0.6Mq for all cases. The dotted line is the L-T c relation of 
an isolated cooling DA WD with pure layers of mass M H = 6 X KT 5 M and 
M Hc = 6 x lfr 3 M Q Mansenl 19991 . 



The dotted line is a rough approximation found by identi- 
fying r max with the core temperature of a cooling WD, giving 
L <x T^ X M, and by assuming that those losses are matched 
by compressional heating, so that L » ?>(M)kT m!lx / /im p (from 
Appendix [E}. Putting these together and using the point 
at (M) = 10" n M Q yr" 1 to set the constant gives L((M)) = 
1.08 x 10- 4 L Q ((M}/10- n M Q yr 1 ) 5 / 3 (Af/0.6M Q r 2 / 3 (dotted 
line). Calibrating the cooling luminosity separately from the 
(M) = curve on the earlier Figure Ogives a prefactor which 
differs from this by less than a factor of 2, demonstrating that 
the estimate from Appendix lElis good for low (M)'s. While 
the L oc (M) 5 ' 3 dependence is fairly well matched at low 
(M)'s, the full calculation gives increasingly lower luminosi- 
ties at higher (M)'s because when fusion provides some of 
the energy, T m . dx need not be as high, so that the luminosity is 
less. The lower solid line gives L(T C — > 0) for M = 1 .0M Q and 
is fairly consistent with (1.0/0.6)" 2 / 3 times the M = 0.6M Q 
curve, shown as the dashed line. The M = 1.0M Q curves fol- 
lows L oc (M) 5 ' 3 to higher (M) than that for M = 0.6M Q . 

Finally we explore the dependence of the L-T c relation on 
M acc . Generally, the energy released by compression in the 
envelope should increase with M acc due to the lower n for 
H/He relative to C/O as it appears in the compressional heat- 
ing term (see eq.[3j. Energy from hydrogen fusion should also 
increase with M acc . The upper curves in Figure[3]shows L for 
M acc of 0.5, 1, 3, and 5 x 1O _4 M for a 0.6M Q WD accreting 
at (M) = 10" 10 - 4 M Q yr -1 . This accretion rate, w 4 x 10" I0 M Q 
yr" 1 , is characteristic of gravitational wave driven accretio n 
for dwarf novae below the period gap ( Kolb & Baraffe 1999). 
Also shown for comparison is (M) = 1O _1I M yr -1 . As M acc 
is increased from 0.5 to 1 x 10~ 4 M Q , three important features 
are apparent. At T c < 3 x 10 6 K the exiting surface luminosity 
L is nearly unchanged because the additional energy released 
in the envelope is transported into the core with a negligible 
rise in r max . For 3 x 10 6 K< T c < 10 7 K, T c is high enough 



that the additional compressional heating energy is forced out- 
ward and acts as an additive factor on top of the cooling lu- 
minosity and there is a small increase across this range. At 
T c > 1.5 x 10 7 K, the density has become high enough that 
the large amount of energy released by fusion has destabilized 
the envelope at these temperatures. This is discussed more in 
Section[3] 

The crossover behavior in T c between core heating at low T c 
and cooling at higher T c can be seen directly in Figure@]which 
shows L core , the luminosity exiting the core into the base of 
the accreted layer. In order to maintain a logarithmic scale, 
two panes are displayed, the upper displays positive L core and 
therefore core cooling, and the lower -L COK to show core heat- 
ing. When increasing M acc from 0.5 to 1 x IO^Mq (solid 
and dashed curves respectively) the core heating increases at 
T c < 4 x 10 6 K. The core cooling, present at higher temper- 
atures, remains fixed with this M acc increase; the additional 
compressional heating is increasing L. Note that at a single 
T c ~ 4-5 x 10 6 K, both core-heating and core-cooling enve- 
lope states are available, depending on the value of M acc . 

Our survey of parameters is now complete, and in Section 
we use L core to discuss the thermal evolution of the C/O 
core. Of the parameters discussed, (M) will be an independent 
variable while T c will be found from (M) and M. However, 
M acc has a maximum value set by the stability of the envelope, 
which we now discuss. 

3. CALCULATING THE IGNITION MASS 

The accreted mass, M acc , increases until a thermonuclear 
instability causes a classical nova explosion at M acc = M; gn , 
which depends principally on the temperature in the accreted 
layer. There are two ignition modes: (1) A traditional thin 
shell instability created by the temperature sensitivity of ther- 
monuclear reactions, or (2) a global instability of the envelope 
relevant when K con d <C Kf-i (r<8xl0 6 K here). The mode 
giving the lower M; gn is the active mode, and each is elabo- 
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FIG. 4. — Luminosity, L core , entering the accreted layer from the core below 
as a function of T c for various accumulated layer masses M aC c on a 0.6Mq 
WD. The top panel shows the rate of core cooling (L C ore) and the bottom the 
rate of heating (— Lcore). The accretion rate is set at (M) = 10~ wa Mq yr -1 , 
and A/acc takes values of 0.5, 1, 3, and 5 X ICT^Mq as indicated. For clarity 
the curves have been truncated at the turning point seen at high T c in Figure 
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rated below. 

The instability manifests in the L-T c relation as a turning 
point, or a maximum T c as a function of L, as seen earlier in 
Figure [3] A general discussion of the relationship between 
turnin g po i nts a nd stability in thin shell accretion appears in 
iPaczvnskl ( 119831) . Due to the large heat capacity of the core 
and the choice of the lower boundary as described in Sec- 
tion EO] T c is essentially constant while M acc increases and 
the turning point moves to lower temperatures. If we de- 
fine M[ gn (T c ) to be the M acc at which the turning point hap- 
pens at T c , the envelope structures with M acc < Mi gn (r e ) form 
a connected quasi-static sequence at constant T c terminating 
at Mi gn . Evolution forward in time from this envelope state 
must be followed in a fully time dependent fashion rather 
than quasi-statically, and so we consider this terminus the ig- 
nition point of the classical nova explosion. The envelope is 
assumed to have switched to a "high", burning state about 
which will make two assumptions: (1) it is short-lived and 
therefore of small thermal importance to the core and (2) dur- 
ing it, most or all of the accreted material is ejected. The thick 
solid curves in Figure|5]show the conditions at the base of the 
accreted layer at ignition for the models equilibrated in the 
manner discussed in 5j4]f° r two values of M. Two values of 
Xi He are shown at each M, representing the range expected in 
cataclysmic variable binaries. 

The lower dot-dashed line in Figure|5]gives the divison be- 
tween radiative and conductive heat transport, where Kf-f = 
^cond- Above this line radiation dominates the heat transfer at 
the base of the accreted layer and the instability is well pre- 
dicted by a traditional thin shell instability criterion. Since the 
scale height of the envelope h = P / pg <C R, it has a positive 
heat capacity and is therefore subject to a thermal instability 
which occurs when the inequality 
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FIG. 5 . — Conditions at the base of the H/He accreted layer at ignition of a 
classical nova for M = 0.6Mq and M = 1 .2Mq . The thin solid curves are the 
local condition given by equation the dotted curves are the global ignition 
condition given by L Q (7i/10 8 K) 2 5 (M/0.6M Q ) = M iga e N (p b ,T b ), and the 
thick curves are the actual conditions found for equilibrated mo dels. The 
dashed lines are the base conditions from ignition calculations bv lFuiimotol 
(1982). 



is satisfied at the base of the accreted material. Here 
e nuc is the energy generated by the fast component of 
nuclear burning; the important chains are l2 C(p,j) l3 N 
(the first step of CNO), 1 H(/7,/3 + ^)D(/?, 7 ) 3 He (the first 
two steps of the pp chains), 3 He( 3 He,2p ) 4 He, and 
3 He(a,7) 7 Be(e-,^) 7 Li(p,a) 4 He llClavtonlll983l . To repre- 
sent local cooling due to heat transfer we define 



^cool 



g 2 4acT 4 
P 2 3k 



(6) 



reflecting the rate at which a thermal perturbation on the size 
of the pressure scale height, h, damps out. The use of divi- 
sion by the scale height in lieu of a derivative to f orm e cnn i is 
comm only termed a one-zone approximation. See Paczvhski 
(1983) for an in-depth example. Equation (Cf) further ap- 
proximates that the pressure response to a small tempera- 
ture perturbation is negligible, a good approximation when 
h -C R. This ignition condit i on is identical to that used by 
Fuiimoto, Hanawa, & Mivaii (1981) to study X-ray bursts, 
though our notation is somewhat different. When the inequal- 
ity of equation (|5} is satisfied, thermal perturbations will grow 
since there is not sufficient cooling to halt the nuclear run- 
away. We evaluate the derivatives using our full equation of 
state and k to obtain a maximum temperature at a given den- 
sity for the base of H/He layer, shown as the thin solid curves 
in Figure [3 for the two values of X) He . When the base tem- 
perature in the H/He layer exceeds this value, ignition occurs. 
This agrees with the location of the turning point in the L-T c 
relations in the radiative regime (See Figure |2J. 

In the radiative regime the opacity is due to free-free and 
at moderate temperatures the triggering nuclear reaction is 
3 He+ 3 He fusion. Both the heating and the cooling can be 
well approximated by power-laws, so that the ignition con- 
dition can be written in closed form. Expanding e nuc about 
T = 10 7 K and using an ideal gas equation of state, equation 
(|5} becomes 
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where T 7 = T/10 7 K and g & = g/lO s cm s , matching the 
radiative portion of curves shown in Figure|5] 

When the base of the envelope becomes degenerate, i.e. 
when K con d <C Kt-f, the turning point happens at a lower M acc 
than the one-zone model predicts, as shown by the divergence 
of the thick lines in Figure [5] from the thin lines of the one- 
zone prediction. The second ignition mode has become im- 
portant and involves an interaction between the radiative and 
conductive regions in the envelope. In order to clarify the en- 
velope state in this global ignition mode, we now present how 
this maximum T c arises from a simple analytical model when 
the base of the accreted layer is in the conductive region. De- 
fine Tb as the temperature at the base of the accreted layer; the 
location of the point of transition to conduction at the bottom 
of the radiative layer leads to L = Lo(7i/7b) 2 ' 5 , where Lq and 
Tq are fiducial parameters. Using a coefficient Kb to character- 
ize the thermal conduction, the heat coming from the core is 
L CO re = Kb(T c - Tb) I zi,c where Zbc is the distance below the base 
of the accreted layer at which T = T c and we ignore the depen- 
dence of K on temperature and density. If the nuclear heating 
during accumulation is approximated by = e^ipb, 7i)M acc , 
where pb is the density at the base of the accreted layer de- 
rived from Pi, = gMi gn /4irR 2 , these quantities are related by 
L = L cois + Ln, neglecting compressional heating, or 
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Here we have taken Ti, = T()(L/Lq) 1 I 2 - 5 and expanded about 
T = 7 x 10 6 K so as to quantify the temperature dependence of 
the burning rate in the regime of interest; only the p+ p — > 3 He 
reaction chain is included at this low temperature. For low 
luminosities, T c increases with L, but eventually the nuclear 
heating becomes strong enough that a maximum T c is reached. 
With the outer temperature boundary of the conductive region 
set by the radiative-conductive transition and the inner tem- 
perature boundary fixed by the thermal inertia of the core, the 
heat released by nuclear processes cannot be transported out 
of the conductive region, and therefore builds up rapidly, lead- 
ing to an explosion. 

A simple approximation of the maximum of equation 
(|8j is obtained by finding when the opposing terms, L 
and M acc ejv, are comparable. This gives the simple form 
io(7l/r ) 2 - 5 (M/0.6M Q ) = M ign e N (p b , T h ). This expresses that 
Mi gn is the accumulated mass for which the outgoing lu- 
minosity is wholly provided by nuclear burning, giving a 
maximal envelope state before CN runaway. By expanding 
£n at the lowest T point of the solid line for M = 0.6M Q , 
T = 4.36 x 10 6 K, p = 5.7 x 10 3 gr cm" 3 , using only the pp 
reaction as above but including the screening correction, we 
obtain oc T 5J3 p l 11 . The dotted line shown in Figure|5]is 
this form with the prefactor adjusted to fit the lowest temper- 
ature point of our actual envelope calculation for M = 0.6M Q 
(the lower thick solid lines). Using an L-Tb relation fitted to 
the (M) = line in Figure[Ogives a result that is parallel and 
only a factor of 2 lower in density, mostly due to the neglect- 
ing the inward-bound heat from the dropped term. 

The dashed l ines in Figure [5] give the results found by 
Fujimoto ( 1982) for M = 0.6M Q (lower line) and 1 .0M Q (up- 
per line). Our results are consistent with these in the demon- 
stration of two ignition modes. The largest difference is due to 
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FIG. 6. — Luminosity exiting the WD core as a function of M acc for several 
values of T c as indicated in units of 10 6 K. All curves have M = 0.6Mq and 
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our inclusion of 3 He in the accreted matter as a nuclear energy 
source. Our calculation of ignition w i th Xs = is roughly 5% 
higher in temperature than Fuiimoto ( 1982) in the radiation- 
transport dominated regime. Much of this 5% discrepancy is 
likely from our updated opacities. 

4. THE THERMAL EQUILIBRIUM OF THE WHITE DWARF CORE 

Section |2] introduced the envelope models which pro- 
vide L(M acc , (M), T C ,M). These same models also provide 
L core (M acc ,(M),T c ,M), the luminosity passing through the 
base of the envelope from the core. The variation of L core 
with M acc at several r c 's and two accretion rates is shown in 
Figure[6] Each of these curves extend in M acc to the CN igni- 
tion 1 and make clear the amount of heat exchanged between 
the envelope and the WD core over the CN cycle, with core 
cooling occurring at low M acc and core heating at high M acc . 
This allows us to calculate an equilibrium core temperature, 
Tc.eq, where the heating and cooling of the core are equal over 
the CN cycle (Townsley & Bildsten 2002). At the onset of ac- 
cretion in the binary, the WD core has T c < T c eq and is heated 
during the CN cycle. We now find the time it takes to reach 
r c , eq at a given (M). 

At low (M) any long-term compression (expansion) of the 
WD core due to mass gain (loss) has little impact on the enve- 
lope (see Appendix A). Hence, for simplicity, we assume that 
the accumulated mass is ejected in a CN, keeping M constant. 
We begin by defining an average core heating rate presuming 
that M and T c are constant while M acc changes, 



,)((M),T C ,M): 



1 

tcN 



^core 

(M acc (t),(M),T c ,M)dt 



(9) 

Using this (L core ) to make statements about the WD's evo- 
lution implicitly assumes that the star's only memory of the 
CN cycle is the net heat gained or lost during it, otherwise it 
returns to an identical state to that at the same stage in the pre- 
vious cycle. Figure shows the dependence of (L C0le ) on T c 

1 The CN explosion has a negligible effect on the thermal content of the 
core due to its short duration compared to the envelope accumulation time. 
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FIG. 7. — Time averaged luminosity escaping from the core during the 
accretion of matter up to a classical nova event. The curves have (M) as 
indicated in units of Mq yr~' and Mwd = O.6M0. 



for several accretion rates. We find that (L core ) increases with 
T c from two effects: 1) The cooling luminosity is greater at 
higher T c , causing stronger core cooling for small M acc , and 2) 
the maximum M acc is larger at lower T c , allowing for a longer 
period of core heating when M acc is large. The equilibrium 
point where (L cole ) = forM = 0.6M Q and (M) = 1(T 10 ' 4 M Q 
yr" 1 , is T c « 5.5 x 10 6 K. For T c rj r ceq , the envelope transi- 
tions from a state which allows the core to cool at low M acc to 
a state which heats the core at high M acc . The positive slope of 
(L core ) at T c w r c eq implies that the equilibrium state is stable, 
i.e. a higher T c at the same (M) causes the core to radiate and 
lower its temperature. The figure also makes clear that r c , eq 
increases with (M) in order to maintain equilibrium under the 
additional heat deposition. 

The specific heat of the liquid WD core is w 3kJ fxim p , 
where /i, w 14 is the ion mean molecular weight (see Ap- 
pendix A). When T c < r c eq , the core is heated by the envelope 
at the rate shown in FigureQ allowing for a simple integration 

<M (10) 
Him p dt 

to yield the time to reach the equilibrium. For a compari- 
son which does not depend on the initial T c , we approximate 
{Lcore} ~ (T c -T ce q)(dL/dT c \) where the derivative is evalu- 
ated at the equilibrium point, giving the form T c - T c eq = (T[ n - 
Tc,w)e~'l T where r = 3kM / fiim p (dL/dT c )o. This e-folding 
time is 0.22, 0.32, 0.45 and 1.9 Gyr for (M) = 10" 8 , 10" 9 , 
10~ 10 , and 1O _I1 M yr" 1 respectively. If (M) changes on a 
timescale shorter than this, T c will no t be maintained at its 
equilibrium value. In previous work iBildsten & TownslevI 
120031) we estimated (dL/dT c ) fa 3k(M)/2^ e m p . This esti- 
mate matches the results from the full envelope calculations at 
(M) ~ fewx 10" n M Q yr" 1 , but the actual (M) dependence is 
much weaker than estimated due to complications introduced 
by the varying ignition mechanism and nuclear energy source. 
For example, heating of the core is relatively less efficient at 
(M) >2x 1O" 1O M yr" 1 when the base of the accreted layer 
never becomes conductive. 



The reheating time at a constant (M) from a reasonable 
starting temperature can be calculated directly from equation 
JlOi and the curves in Figure for all but the highest accre- 
tion rates. The time for the core of a 0.6M Q C/O WD to 
heat u p from 2.5 x 10 6 K , corresponding to a cooling age 
4 Gyr (Salaris et al. 2000), to 90% of the equilibrium tem- 
perature is 0.7 Gyr, 0.9 Gyr and 2.4 Gyr for 10" 9 , 10" 10 and 
10" 1 1 Mq yr" 1 respectively. The mass that must be transferred 
in the binary in each of these cases is thus 0.7, 0.09, and 
0.024M Q . The latter values are quite reasonable expected 
fluxes in terms of evolution of cataclysmic variables (see 
Howell, Nelson, & Rappaport 2001), and thus WDs in CVs 
with 1.3 < P 01 b < 2 hours are expected to be in equilibrium. 
However, the flux necessary at 10~ 9 Mq yr" 1 is quite high, so 
that systems with P or b > 3 hours are unlikely to have reached 
equilibrium. This will be mitigated to some degree by the 
fact that all of the ~ 0.4M Q transferred above the period gap 
occur s at (M) > 1O" 9 M yr" 1 jHowell. Nelson. & Rappaport 
1200 11) . which has stronger core heating. A system which 
comes into contact at P or b « 6 hours should be coming ap- 
proximately into equilibrium when P or b ~ 3 hours, just above 
the period gap. We defer a detailed calculation and discussion 
relating the evolution of the interacting binary and the t hermal 
state of the WD to a companion paper ( Tow nslev & Bildstenl 
l200H) . 

5. CLASSICAL NOVAE IGNITION MASSES 

The principal outcome of our calculations is the equilibrium 
temperature at which the heating and cooling of the WD core 
are balanced over the CN cycle. Once this equilibrium has 
been reached, observable quantities do not depend on the age 
of the WD, but rather on M and (M), which are determined 
by the evolution of the binary system. 

The top panel of Figure [8] shows T c , eq for a range of 
(M). Curves for three values of M are shown, two for 
each with Xi He = 0.005 and 0.001, the upper and lower 
limits on expected values for the accreted mate rial (e.g. 
iD'Antona & Mazzitellill98lllben & Tutukovll984 . The de- 
pendence of T c eq on these parameters is complicated due to 
the fact that two transitions are taking place over this broad 
range of accretion rates. There are three regimes: (1) At 
(M) > 3 x 10~ 10 M Q yr" 1 , the thermonuclear runaway is well 
described by a one-zone model, and the dominant nuclear 
fuel triggering the explosion is 3 He, T c , eq increases with (M) 
in order to balance the higher rate of heat deposition in the 
envelope. Larger Xi He causes the CN explosion to happen 
sooner, thus truncating the heating phase relative to lower 
Z3 He and leading to a lower r c eq . The lower limit of this 
regime depends on M since higher gravity provides a more 
stable envelope in the one-zone model. (2) At 3 x 10" 11 < 
{M)/(Mq yr" 1 ) < 3 x 10" ll) the one-zone model of ignition is 
no longer sufficient (see Section[5}, but the dominant nuclear 
heat source is still 3 He. The contrast in T c , eq with M has be- 
come quite small, while the dependence on Xs He remains. (3) 
At (M) < 3 x 10" n M Q yr" 1 , the dependence on X? He is mostly 
gone and the dominant nuclear heating source is pp. 

While T c is not directly observable, the resulting ignition 
masses, M; gn , can be compared to measured ejection masses. 
The bottom panel of Figure|8]shows M; gn for the equilibrated 
stars. This is the first self-consistent evaluation of the CN ig- 
nition mass since we use our calculated T c , ecj ((M}). The igni- 
tion mass decreases as the mass of the WD increases. This 
dependence is familiar from previous work, but is usually 
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FIG. 8. — Equlibrium core temperature and mass of the accreted layer at 
classical nova ignition, M\ gn as a function of time averaged accretion rate, 

(M). Curves are shown for three WD masses, M, and two different mass 
fractions of 3 He at each M. 



stated as the contrast expected for a constant ignition pres- 
sure, P ign w GMMi gn /4ttR 4 , yielding M iga oc R 4 /M. Clearly 
the ignition pressure for a fixed mass is not constant, as T c 
changes with (M), thus the simple scaling does not hold. 

The contrast in nuclear energy contributions is clearly 
demonstrated; 3 He triggering has a pronounced effect on M; en 
at (M) > 1O- 1O M yr" 1 and barely any for (M) < 4 x 1O~ U M 
yr -1 where pp is the dominant contributer. Our M lm values 
with X 3 HP = 0.001 are 1/2 to 1/4 of those found by Fuiimoto 
( 1982) (see his Figure 7), the fraction decreasing with in- 
creasing (M). This difference is surprisingly small given that 
his (M)'s were order of magnitude estimates, in contrast to 
our self-consistent calculatio ns. Th e grid of CN simulation 
results in Prial nik & Kovetzj i 19951) includes values of M- lgn 
(their w acc ) for M = 0.65, 1 .0 and 1 .25M Q , (M) = 10" 10 , 10" 9 
and 10 _8 M o yr" 1 with T c = 10 7 K. For M = 0.6 and 1.0M Q 
at 10" 8 and 1O _9 M yr" 1 , we find about 1/2 of their M ign for 
Xi He = 0.001. Whereas at 10- 10 M Q yr" 1 , again with X, ne = 
0.001, we find 1, 1.4, and 3.1 times their M ign atM = 0.6, 1.0, 
and 1 .2Mq respectively. We understand these differences as 
follows. At (M) = 10- 10 M Q yr 1 our equilibrium temperature 
is less than 10 7 K thus leading to larger ignition masses. Sec- 
ond, our inclusion of 3 He in the accreted material decreases 
Mi gn with respect to calculations which include only CNO en- 
richment leading to our lower M lgn values at the higher accre- 
tion rates where T c eq is close to 10 7 K, and somewhat coun- 
teracting the affect of the lower T c at (M) = 1O" 1O M yr" 1 . 

Our calculations of M lgn agree quite well with those of 
iMacDonaldl (119841) for (M) < 10- 105 M o yr" 1 . This is where 
pp burning dominates and the base of the accreted layer is 
degenerate at ignition. However, his r c e q's were not reported. 
At larger (M) ~ 10- 9 -10 _8 M Q yr 1 our M; gn are nearly a factor 
of 10 lower than MacDonald. This discrepancy is larger than 
can be explained by 3 He ignition. A s discusse d in MacDonald 
( 119831) . the "cold" models utilized in lMacDonaldl(ll984l) leave 
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FIG. 9. — Comparison of our predicted classical nova ignition masses 
(Mj gn ) to measurements of classical nova ejected masses for systems which 
have a measured orbital period in Warner (2002). The ejected mass references 
are: V1974 Cyg (Gehrz et al. 1998), QU Vul (Shin et al. 1998), V1668 Cyg 
(Gehrz et al. 1998), V1425 Aql (Lyke et al. 2001), V838 Her (Gehrz et 
al. 1998). The two lines at each M as indicated are for Xt He = 0.001 and 
0.005, giving larger and smaller M„ n respec tively. At low peri ods (M) is that 
expected from gravitational radiation I Kolb & Baraffe 1999) and at P or i, = 
3 and 6 hours we use (M) = 10~ 9 and 10~ 8 Mq yr~' respectively, which is 
approximately that expecte d for above the period gap due to magnetic braking 
IHowell. Nelson. & Rappaportl200il) . 



out one of three ignition conditions. The omitted condition 
turns out to be quite similar to the one-zone condition and 
therefore essential at these higher (M)'s. 

Accumulated masses are not directly measurable, but the 
mass ejected in a CN, M e j, is known for many systems. 
A robust comparison between Mj gn and M ej would tell us 
whether the WD mass is increasing or decreasing on long 
timescales during accretion. Because Mj g „ depends on (M), 
we must estimate (M) for a given CN system, requiring us 
to know the orbital period and to make some presumption 
about its relationship to (M). For P or b < 2 hours, (M) is that 
give n by angular mo mentum loss due to gravitational radia- 
tion llKolb & Baraffel U999l) . For f or b > 3 hours mass trans- 
fer is driven by magnetic brak ing. To approximate the stan- 
dard model as presented in e.g. lHowell. Nelson. & Rappaport 
( 1200 ll) . we have used two values, (M) = 1O" 9 M yr" 1 at 
P orb = 3 hours and (M) = 10~ 8 M Q yr -1 at f rb = 6 hours. A 
cross-comparison between M; gn and M e j is shown in Figure 
|9l which shows M ign for M = 0.6 and 1.0M Q and Xi He = 
0.001 and 0.005. The five M e j points, which have large er- 
rors, are those systems where both the orbital period and 
ejected mass are known. One of these systems, VI 974 
Cyg, has a kinematic WD mass estimate of 0.75- 1.07M Q 
( iRetter. Leibowitz. & Ofeklll997l) . Other estimates are less 
certain. Due to the small number of points, large errors, and 
lack of measured WD masses, the comparison to the current 
data-set is not definitive, but does roughly agree, pointing out 
that most of the accreted mass is ejected during the CN. We 
have not found overwhelming evidence for excavation of the 
WD. 

There are now over 50 CN with measured orbital periods 
(Warner 2002). Our ignition masses will allow for a calcu- 



10 




Accretion Rate (M /yr) 



FIG. 10. — Ranges of WD T e g as a function of accretion rate for WD 
mass M = 0.6, 1.0, and 1.2Mq (bottom to top). The minimum is evaluated 
when M acc = 0.05Mi en and the maximum when M ;lcc = 0.95Mj gn . For com- 
pa rison, in 1 Gyr a 0.6 (1. 0)Mp, isolated WD cools to T cff = 9 X 10 3 (14 X 10 3 ) 
K ISalaris et all200(l) . 



lation of the CN orbital period distribution when combined 
with CV evolution. For example, the tenfold contrast in ig- 
nition mass above and below the period gap makes the CN 
occurence rate nearly a factor of 100 lower below the period 
gap than above. Currently 10% of the CN are below the pe- 
riod gap, if these systems all suffer the same selection effects 
then the implied population below the period gap is 10 times 
the number of systems above the period gap. 

Our envelope models also predict r et f = (L/4-TTR 2 a SB y/ 4 
which increases with M acc . Figure [lOl shows, as a function 
of (M) for T c = r e eq , the range of T e ff traversed between 
M acc = 0.05M ign and 0.95M ign for M = 0.6, 1.0, and 1.2M Q . 
We find that the surface flux L does not depend strongly on 
mass, and thus roughly T e g oc R~ 1 ' 2 . Short horizontal lines 
with this dependence are shown for comparison. The WD 
r e ff has been measured for a number of Dwar f Novae systems 
when the disk is in quiescence (ISionlll999l) . and a detailed 
compa rison to such measurement s is being published sepa- 
rately (i^wnsie^&Biidite^lQQl). 

6. CONCLUSION AND DISCUSSION 

A quasi-static WD envelope model has been used to find the 
equilibrium core temperature, T cm of an accreting WD under 
the approximation that exactly the accreted layer is ejected 
in each classical nova outburst. Such an equilibrium can be 
reached in < 1 Gyr, allowing for CVs to come close to equi- 
librium by the time they reach P or b ~ 3 hours, and insuring 
equilibration for P or b < 2 hours. Having solved for r ceq it is 
then possible to predict the CN ignition mass M[ gn and the 
thermal contribution to the WD T e ff. This is the first deter- 
mination of M- isn which utilizes a WD state determined from 
(M) and the cooling and heating of the core during the CN 
buildup, making it fully self consistent. When compared 



with several papers from the CN literature (Fuiir notollT98l 
iPrialnik & Kovet3[l99l lMacDonaldlU98l we find that the 
inclusion of 3 He leads to lower M ign for (M) > 10" I0 M Q yr" 1 , 
and that for (M) below this the particular author's assumption 
concerning T c , which we calculate consistently, is a determin- 
ing factor. When compared with observed ejected masses, our 
Mig n are compatible with M e j ~ M; gn . Had this comparison fa- 
vored either M e j > M; gn or the opposite, the equilibrium would 
take place at some T c such that L core has the value appropri- 
ate for the implied secular change in core mass, which can be 
calculated using the formalism presented in AppendixlAl 

Over twenty Dwarf Novae have been observed in quies- 
cence, when the accretion rate is low and the WD pho- 
tosphere is detected and T e s measured. The theoretical 
work p resente d here is c ompared to these observations in 
iTownslev & Bildstenl (120031) . allowing us to constrain the WD 
mass and the time averaged accretion rate, (M). We show 
there, for systems with P 01 -b < 2 hours, that if (M) is given 
by gravitational radiation losses alone, then the WD masses 
are > 0.8M Q . An alternative conclusion is that the masses are 
closer to 0.6M Q and (M) is 3-4 times larger than that expected 
from gravitational radiation losses. 

It is well known that an isolated WD will pulsate 
when its r e g is in t he approximate range 11000-12000 K 
(Berg eron et alll995l) . While a difference in the atmospheric 
composition, H/He mixture under accretion versus pure hy- 
drogen on the non-accreting case, will shift this range, it is 
likely that a similar pulsation mechanism will be active in ac- 
creting WDs. Our calculations indicate that accreting WDs 
with M = 0.6-1.0M Q should be near this range when (M) = 
fewxlO" u M yr -1 (See Figure [TuTi. This (M) is typical 
of that expected in CVs when accretion is driven by emis- 
sion o f gravitational radiation, P or b < 2 hours ( iKolb & Baraffd 
1999). In fact one system, GW Lib, has been found which 
does exhibit precisely thi s type of variability ( Szko dv" ;t alJ 
120021 Ivan Zvl et afll2000l) . Using the interior models devel- 
oped here, we are now undertaking a seismological study of 
these systems. This offers the tantalizing possibility of facili- 
tating measurement of the WD mass and spin and determina- 
tion of gross internal structural features such as the size of the 
accreted layer. 

We thank Ed Sion for numerous conversations, Paula 
Szkody for up to date information on the DN observations, 
Brad Hansen for insights on modern WD modeling, Francesca 
D'Antona for questioning our equilibria, and the referee for 
constructive criticism. This research was supported by the 
National Science Foundation under Grants PHY99-07949 and 
AST02-05956. Support for this work was provided by NASA 
through grant AR-09517.01-A from STScI, which is operated 
by AURA, Inc., under NASA contract NAS5-26555. D. T. is 
an NSF Graduate Fellow and L. B. is a Cottrell Scholar of the 
Research Corporation. 



APPENDIX 

A. THE THERMAL STATE OF THE UNDERLYING WHITE DWARF 

The thermal state of the WD core can affect the envelope structure through the heat flow into or out of it on the evolutionary 
times we are considering. Due to varying ejected mass in classical novae outbursts, we do not actually know whether the WD 
core is increasing or decreasing in mass. However, we show here that for a C/O core, the compression or decompression rate of 
the WD core will have little impact on the overall envelope calculations. This is less clear for tho se WDs with pure helium cores 
that are likely pr esent in binaries below the pe r iod gap (e.g.lHowell. Nelson 1 & Rappapo rtl200ll) . Our discussion follows closely 
the insights from Nomoto & Sugi motol ( fl977i) . lNomotol(ll982l) . and lHernanz et alJ ( 119881) 
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The WDs of interest here consist mostly of 12 C and ls O so we construct the deep interior (far beneath the freshly accreted 
hydrogen and helium) using the degenerate electron equation of state for \i e = 2, neglecting the small 22 Ne fraction. The ions 
completely dominate the thermodynamics and thermal conductivity of the liquid WD interior. For a classical one component 
plasma (OCP) with ion separation, a, defined by a 3 = 3/47rn,, where n, = p/Am p , the importance of Coulomb physics for the ions 
is measured by 

(Zef 1/3 /10 7 K\ /Z\ 2 /16\ I/3 

^r=57.7p 6 —Jig) T 



where p& = p/10 6 grcm 3 . We will work in the liquid regime, V < 175 ( Potek hin & Chabrierll2000l and references therein). 
For a 0.6(1. 2)M Q pure Oxygen WD, this requires that the core temperature, T c , exceed s 5(16) x 10 6 K (this is within the 
range expected for the low mass WDs from section 5). Following the recent work of iDeWitt. Slatterv. & Chabrierl i 19961) and 
Chab rier & Potekhir] dl998l) . we used the linear sum rule to evaluate the Coulomb energy in the mixed liquid phase and the 
resulting specific heat, cy. We then integrated throughout the stellar model to fi nd the total heat capacity of the white dwarf, 
Cv = J cydM, as a function of the isothermal WD temperature. In agreement with Potekhin & Chabrier (2000), we find that Cy 
is nearly that of a classical crystal, Cy = 3k{M / piin p ), for temperatures in the range of the melting value and slightly above and 
below. 

Elec tron conduction dominates the heat transport in the deep liquid WD interior and we use the work of lYakovlev & Urpir] 
( 1980) for our estimates of the thermal conductivity, K, when electron-ion collisions dominate. For an isothermal WD, the time 
it takes for conduction to transport heat from the outside edge, R, to an interior radius r is 

x 2 




^ = 4^ \jk) pdr ) < (A2) 

as shown by Henvev & L'Ecuver ( 1969). This is just a slightly more sophisticated version of what Nomoto & Sugimoto fl97l 
call Th = cypHp/K, the time for heat transport across a pressure scale height, H p . Since cy is nearly independent of temperature, 
the main T scaling of f conc | is from the conductivity, yielding f con d ocl/T. Figure lATTl plots f con( j as a function of the mass location 
within pure, isothermal (T = 10 7 K) carbon WDs of masses M = 0.4,0.6,0.8, 1.0 and 1.2 M Q . The timescale is only slightly 
different for pure oxygen WD's. 

The compression of the WD interior is adiabatic if accretion compresses it on a timescale much shorter than f conc j ss (2-4) x 
10 7 (1 7 K/r) yr. To relate the compression rate to M, we start at the WD core and write dlnp c /dt = {din p c / dlnM)dlnM / dt 
(e.g. lNomotoi ri982). For a low-mass WD, a = dlnp c /dlnM = 2, but a increases as M approaches the Chandrasekhar mass, 
having values of a = 3.27,4.26,6.17, and 10.71 for M = 0.6,0.8, 1.0, and 1.2M Q . For example, the accretion rate required to 
adiabatically compress the interior of a 0.8M Q WD when T = 10 7 K is > 6 x 10~ 9 M Q yr" 1 . Simulations for rates this rapid have 
been carried out (e.g. Nomoto & Iben 1985) and clearly show the central adiabatic compressio n. Since the ions do minate the 
entropy there, we calculate the adiabat in the liquid state using the internal energy given by Chabrier & Potekhh] (EH) to find 

(d\nT\ 0.91+0.14r 1/3 

\dh^) ad ~ 1.22 + 0.4ir 1/3 ' ( } 

within 1% of the adiabat given bv lHernanz et al.l dl988l) for T > 1 to crystallization. For accretion rates where the compression 
time and conduction time are comparable, there is clearly time for heat to be transported, but not adequate time for maintaining 
an isothermal core. 

However, for M<6x 1O~ 9 M yr" 1 , there is time for the core to remain isothermal under the compression (Hernanz et al. 
1988), simplifying the treatment of the WD core for the M's of most interest to us. Unlike our treatment of the atmosphere, for 
the deep interior, we need to account for the changing global structure of the WD. Nomoto (1982) showed that this is greatly 
simplified by rewriting the entropy equation in terms of the variable q = M r /M, where M r is the mass interior to r and M is the 
total mass. Under accretion, we relate q/q = M r /M r -M/M, and rewrite the time derivative as d/dt = d/dt \ q + qd/dq\ t . In the 
entropy equation, we want the Lagrangian time derivative djdt\u r , so that q = —qM/M, and the entropy equation simply becomes 

ds\ „d\nM f ds \ dL r 



T \dt) q T dt \dlnq), dM r ' (A4) 
under the action of compression. Integrating this equation tells us the luminosity exiting the surface of the WD core 

dlnM f M ( ds \ f M fds\ 

as clearly written by Nomoto & Sugimoto (1977). 
The advantage to writing the terms in this manner becomes apparent once we use 

k 

ds = [adlnT- bd\np\ , (A6) 

fiitn p 

where a = 1.22 + 0.41F 1 / 3 , and b = 0.91 +0.14T 1 / 3 are presumed to be roughly constant throughout the WD . If the ions were an 
ideal gas, a = 3/2 and b=\. Presuming that the WD core is always isothermal at T = T c , then equation jA5> becomes 

u,m„ dlnM f M fdlnp\ d\nT c f M ( dlnp \ dlnM 

L^-^= — / b dM r -aM—^+ / b — — y - ——dM r . (A7) 

kT c dt Jo \dlnq) dt J \dlnM J dt 
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FIG . All. — The thermal conduction time (f CO nd in equation IA2I ) from the exterior of a pure carbon WD to an interior mass point. The curves are for isothermal 
WDs (T = 10 7 K) with masses M = 0.4,0.6,0.8, 1.0 and 1.2M©. 



We rewrite this more clearly by regrouping terms to obtain 

aMk dT r bkT r d\nM 



L+- 



Hirrip dt 



liim p 



dt 



M 



( dlnp\ ( dlnp 
\d\nq) l \d\nM 



dM r . 



(A8) 



The term c/lnp/c/lnM^ is the generalization of the piece discussed earlier in the context of the central point. It gives the 
compression rate at a fixed mass coordinate, q, due to accretion. No moto (1982) showed that it is independent of q within the 
deep interior, so that a = dlnp c /dlnM w d\n p / d\nM\ q , and equation (IA8> becomes 



L+- 



aMk dT r abkT r dM bkT r dhxM 



\iim p dt Pitn p dt p.i>n p dt 



M r d\np, 



(A9) 



where the last term gets large near the surface and is evaluated in Appendix B as part of the outer envelope integration (see 
equation B3). At V w 80, a = 3 and b = 1.5, so for a M = 0.8M Q WD we get the simple relation for compression of the deep 
interior of the WD 

aMk dT c 6AkT c dM 

L+ -£ = e - — , (A10) 

fiiiiip dt Pitrip dt 

implying that the core can remain at a fixed temperature during compression from accretion as long as the exiting luminosity is 

6AkT c dM 



fAitrip dt 



(All) 
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Comparison of this with the corresponding estimate for the accreted layers from AppendixFBl equation iBAl . demonstrates that, 
even if the WD core is accumulating mass on a secular timescale (i.e. no ejection of matter), then, since fit « 14, such an exiting 
luminosity as L C0le would be a small perturbation on the thermal state of the envelope. This statement would change for a helium 
WD, where /i, = 4, in which case the core luminosity could modify the envelope structure. 

B. ANALYTICAL ENVELOPE INTEGRATION 

Here we present the integration of equation Q from Section fusing a simple model for a WD consisting of a radiative outer 
envelope and a degenerate deeper layer that is isothermal. Dropping the nuclear burning and the time derivative term from 
equation leads to the integral 

L={M)J o T-dP. (Bl) 

To allow a simple integration, in the envelope we make the approximation that the luminosity is constant with depth and that the 
opacity has the form k oc pT~ 35 so that P 2 oc T 8 5 . The entropy of an ideal gas is s = k\n(T 3 ' 2 / p)j pm p , which can be expressed 
in terms of pressure alone, allowing evaluation of ds/dP, and 

k ■ C dP kT 
Lem = 0.41 — — (M) / T^- = \.15{M)^- . (B2) 
pm p J P pm p 

Where P 2 oc T has been used to evaluate the integral, the term from the outer boundary has been dropped, and the pressure at 
the inner boundary (the transition to degeneracy) has been written in terms of T c . 

The degenerate section is nearly isothermal and the entropy is completely in the ions, so that s = k\n(T 3 / 2 /p)/pjin p , and the 
density is related to the pressure by P oc p 5 ^ 3 from the surrounding non-relativistic degenerate electrons. From this ds/dP can 
be evaluated, but some limits must be chosen. The lower limit is the pressure at which the transition to degeneracy occurs, 

5 12 

this is simply P tI = T c ' in cgs units for a solar mixture. The upper limit is the pressure at the base of the accreted envelope 
Aase = gM-^/AuR 2 . From our results, at (M) = 10" I0 M Q yr" 1 , T c = 8 x 10 6 K and M ign = 3 x 1O" 4 M . Using these values gives 
P tT = 1 .8 x 10 17 erg cm" 3 and P\, ase = 6.5 x 10 18 erg cm" 3 and thus 

IJ^tM) ( d\nP = 2.\5(M)-^- 

5 pimp J pitn p 



L deg = -—^-(M)ldlnP = 2.l5(M) T7 ^ r . (B3) 



Although the limits taken actually depend on T c , they appear in a logarithm so only the strongest temperature dependence is 
shown explicitly. Using pj pi = 0.6/1.3 for solar material to eliminate p, and adding the luminosity from the nondegenerate and 
degenerate portions of the accreted layer gives 

. 3(M)kT c 



pm p 

This fiducial estimate of the compressional heating proves helpful in our analytic understanding of our initial results. 



(B4) 
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